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The aim of the paper is to present the performance characteristics of a Stirling refrigeration 
cycle in micro/nano scale, in which the working substance of cycle is an ideal Maxwellian 
gas. Due to the quantum boundary effect on the gas particles confined in the finite 
domain, the cycle no longer possesses the condition of perfect regeneration. The inherent 
regenerative losses, the refrigeration heat and coefficient of performance (COP) of the 
cycle are derived. It is found that, for the micro/nano scaled Stirling refrigeration cycle 
devices, the refrigeration heat and COP of cycle all depend on the surface area of the 
system (boundary of cycle) besides the temperature of the heat reservoirs, the volume 
of system and other parameters, while for the macro scaled refrigeration cycle devices, 
the refrigeration heat and COP of cycle are independent of the surface area of the system. 
Variations of the refrigeration heat ratio r R and the COP ratio r £ with the temperature ratio 
r and volume ratio r v for the different surface area ratio r A are examined, which reveals 
the influence of the boundary of cycle on the performance of a micro/nano scaled Stirling 
refrigeration cycle. The results are useful for designing of a micro/nano scaled Stirling cycle 
device and may conduce to confirming experimentally the quantum boundary effect in the 
micro/nano scaled devices. 

© 2008 Elsevier B.V. All rights reserved. 


1. Introduction 

In classical thermodynamics, the cycle models of various heat engines or refrigerators using gases as working substance 
are based on the classical ideal gas equation of state and the classical law for heat transfer between the reservoirs and the 
working substance. For a number of real cycle devices, the predictions of these models about the performance of cycles 
are good enough since in general on the one hand, the temperature inside the working substance is high and its density is 
low enough, and on the other hand, the size of the considered system is large enough. However, at some special conditions 
i.e. sufficiently low temperature or high density conditions, quantum degeneracy of gas particles becomes important and 
causes differentiation of behaviors of the gas cycles since the states of gases deviate from the classical ideal gas equation of 
state. In particular, the system, called quantum gas system, should obey the laws of quantum statistical mechanics instead 
of the classical one. By considering the quantum degeneracy of ideal gases, the performance characteristics of the power or 
refrigeration cycles working with ideal quantum gases have been analyzed [1-6]. Moreover, the finite-time thermodynamic 
theories have been used for analyzing the performance of the quantum heat engines [7-12]. 

Most of the studies about the performance of the gas cycle always neglect the influence of boundary of cycle. That 
is, the performance of the cycle is independent of the boundary of cycle. For a macro scaled gas cycle system, the re¬ 
sults are valid. However, recent advances in nanotechnology [13-15] make the thermodynamic behavior of gases in finite 
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scale a current issue [16-21]. In the literature, it has been found that thermal wavelength of gas atoms confined in a mi¬ 
cro/nano scaled domain may be comparable with the size of the system and the quantum boundary effect of gas particles be¬ 
comes important [16-18]. Especially, based on the boundary effect, a micro/nano thermosize gas cycle, which is the analog of 
thermoelectric generator based on the thermoelectric effect, is designed theoretically and may conduce to confirming ex¬ 
perimentally the existence of the boundary effect [22,23]. Thus, even for the classical gas cycle systems in finite scale, similar 
to the influence of quantum degeneracy on the classical cycle, quantum boundary effect of gas particles in a finite domain 
cannot be neglected and may cause differentiation of behaviors of the gas cycles. And the difference may be measured ex¬ 
perimentally by designing a gas cycles i.e. the Stirling refrigeration cycle devices in micro/nano scale. Then, it is necessary 
to analyze further the performance characteristics of the classical Stirling refrigeration cycle in micro/nano scale, in which 
the quantum boundary effect has not been previously considered. In the present paper, a micro/nano scaled Stirling refrig¬ 
eration cycle working with ideal Maxwellian gases is considered. The inherent regenerative losses, refrigeration heat and 
coefficient of performance of the cycle are derived using the thermodynamic relations of the finite domain. It is seen that 
they are different from those for the general Stirling cycle, which works at a macro scaled system. 

The main purpose of the paper is to present the performance characteristics of a micro/nano scaled Stirling refrigeration 
cycle and analyze the influence of the boundary i.e. the surface area of cycle on the performance of the cycle by some special 
examples. The results attained here provide some guidance for designing of a micro/nano scaled heat exchange device and 
may also conduce to confirming experimentally the quantum boundary effect in the micro/nano scaled devices. 


2. Thermodynamic properties of an ideal gas confined in a finite domain 


Based on Maxwell-Boltzmann statistics, the free energy of an ideal Maxwellian gas confined in a finite domain which 
shape can be rectangle, cylinder or sphere is given [16-18] 
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where L c (I) = h/y8mkt,T is one half of the most probable de Broglie’s wave length of the particles at temperature T, h 
is Planck’s constant, m is the atomic mass and k b is Boltzmann’s constant, N is the total number of particles in the whole 
confinement volume, A is the surface area of the domain and V is the volume of the domain. The first term in Eq. (1) is the 
classical bulk free energy of gas confined in the finite domain, while the second term is the quantum surface energy, which 
unavoidably appears and depends strongly on the boundary i.e. surface area of the domain A when the size of the domain is 
comparable with the mean de Broglie’s wavelength of the particles. Further, these additional control variables on the state 
functions may make the ideal gas confined in a finite domain deviate from its classical behavior i.e. the gas confined in a 
finite domain cannot fill all the space of the domain and the density goes to zero within the layer approached the boundary 
of domain [18]. It is interesting that the nature of gases confined in a finite domain will render a realistic thing to design 
a mechanical gas cycle devices without any physical contact between the piston and cylinder of the cycle. Further, the 
mechanical friction between the piston and cylinder can be decreased largely. In addition, because of the boundary effect 
on the thermodynamic state functions of a gas confined in a finite domain, the performance characteristic of the cycle i.e. a 
Stirling refrigeration cycle in micro/nano scale is different from the macro scaled Stirling refrigeration cycle, such as, the 
changing of the refrigeration heat and the COP of cycle. That is, it may be observed that the quantum boundary effect of gas 
particles in a finite domain will cause the differentiation of behaviors of the classical gas cycles. 

The internal energy and entropy of an ideal gas confined in a finite domain are, respectively, given by 
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where C = ( 2jrmk b ) 3/2 /h 3 . From Eq. (2), one obtains the heat capacity at constant volume as 
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Using Eqs. (3) and (4), one can calculate the amounts of heat exchange in an isothermal and isochoric process. They are given 

by 
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Fig. 1. Temperature-entropy diagram of a micro/nano scaled Stirling refrigeration cycle. 


and 

V f b 3 Nk b A / <— r~\ 

Qat, = J C V (T,A, V ) dl = -Nk b (T b - T a ) + ^3 - (VTb - A,) , ( 6 ) 

where the subscripts a and b refer to the initial and final states, respectively. With the help of the thermodynamic relations 
in the finite domain, we can calculate the coefficient of performance (COP) of a micro/nano scaled Stirling refrigeration 
cycle working with ideal gases and discuss the influences of boundary on the inherent regenerative losses in two isochoric 
processes, the refrigeration heat and the COP of the cycle. 


3. A micro/nano scaled Stirling refrigeration cycle 


A micro/nano scaled Stirling cycle is composed of two isothermal and two isochoric processes. The temperature-entropy 
diagram of the cycle is shown in Fig. 1 , where Q.i 2 and Q34 are the amounts of heat exchange between the working substance 
and the heat reservoirs at temperature T H and T L during the two isothermal processes, Q23 and Q 41 are the amounts of 
heat exchange between the working substance and the regenerator during the two isochoric processes, and V H and V L are, 
respectively, the maximum and minimum volumes and A H and A L are, respectively, the maximum and minimum surface 
areas of the gas system. All heats Q.i 2 , Q34, Q23 and Q 41 are positive. In order to improve the performance of the cycle, a 
regenerator is often used in two isochoric processes. By using Eqs. (5) and ( 6 ), the amounts of heat exchange in the various 
processes of the micro/nano scaled Stirling cycle may be expressed as 
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respectively. From Eqs. (7)-(10), we can obtain the work input per cycle as 
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where r = T L /T H , r v = V H /V L and r A = Ah/A l . Because the net input work must be positive, we have r A < 
r v [1 + 2 yWV L (1 + J^) In r v /L c (T H )A L ]. 

Since the influence of surface area the amount of heat Q 23 flowing into the regenerator in one isochoric process is different 
from the amount of heat Q 41 flowing from the regenerator in the other isochoric process. Furthermore, by using the Eqs. (9) 
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and ( 10 ) the regenerative losses A Qj is given by 


AQj = 023 - 041 = Nk b T H L (S^)-^ (l - -) (1 - Vr) 

4Vjr V L \ r v J y 


( 12 ) 


From Eq. (12) we can observe that the value A Qj is dependent of the volume ratio and the surface area ratio of two isochoric 
processes and can be positive, negative or zero according to the other parameters. 


3.1. Case of AQj > 0 


In the case, according to the Eq. (12), A Qj > 0 means that the range of surface area ratio is r A < r v . 
However, considering the definition of r A , then the range of surface area ratio is only 1 < r A < r v . Further, the amount 
of heat exchange Q 23 flowing into the regenerator in one isochoric process is larger than or equals that of heat exchange 
Q .41 flowing from the regenerator in the other isochoric process. The redundant heat in the regenerator per cycle must be 
released to the cold reservoir in a timely manner because the present Stirling cycle is only operated between the two heat 
reservoirs T H and T L . If not, the temperature of the regenerator would be increased such that the regenerator would not 
operate normally between the T H and T L . Thus, the refrigeration heat Rqb absorbed from the cold reservoir is obtained as 


Rqb — Q34 ~ AQ. — Nk b T H 
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For an ideal gas confined in a macro scaled domain, the second term in Eq. (13) equals zero because L c (T) < < V /A and the 
quantum boundary effect can be neglected. Therefore, the refrigeration heat of classical cycle can be written as 


R c = Nk b T H x lnry. (14) 

However, in micro/nano scale the thermal wavelength L c (T) of gas atoms may be comparable with the size V/A of the 
domain and the quantum boundary effect cannot be neglected. Especially, the quantum boundary effect of gas particles is 
more important as the size of the finite domain is smaller. Therefore, in these cases the second term must be considered. 
From Eqs. (11) and (13), the coefficient of performance of cycle in the case can be expressed as 
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Furthermore, when r A = r Vf the COP of the cycle equals the ideal Carnot COP 


(15) 


Sqb = r /1 - t = e c . (16) 

From Eqs. (13) and (15) we can find that the refrigeration heat and COP of micro/nano scaled cycle are not only dependent 
of the volume ratio of two isochoric processes but also the surface area ratio of two isochoric processes when 1 < r A < r v . 


3.2. Case of AQj < 0 


In the case, similarly, according to the Eq. (12), AQj < 0 means that the range of surface area ratio is r v < 
r A . However, considering the condition W > 0, then the range of surface area ratio is only r v < r A < 

r v [l + 2^/ttVl (l + y/r) In r v /L c (T H )A L ]. Especially, in the case the amount of heat exchange Q23 flowing into the 
regenerator in one isochoric process is smaller than that of heat exchange Q 41 flowing from the regenerator in the other 
isochoric process. Moreover, the inadequate heat in the regenerator per cycle must be compensated from the hot reservoir 
in a timely manner because the present Stirling cycle is only operated between the two heat reservoirs T H and T L , while the 
refrigeration heat Rqb is unvarying when AQ r < 0. If not, the temperature of the regenerator would be changed such that 
the regenerator would not operate normally between the T H and T L . Thus, the refrigeration heat Rq B absorbed from the cold 
reservoir is obtained as 
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From Eqs. (11) and (17), the COP of the cycle can be expressed as 
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In the case, according to the Eqs. (17) and (18), the refrigeration heat and COP of the cycle also depend strongly on the surface 
area ratio of two isochoric processes besides that the volume ratio of cycle and other parameters. Eqs. (13), (15), (17) and 
(18) are the main results of this paper. Further, by using these equations, the performance characteristics of the micro/nano 
scaled Stirling refrigeration cycle may be analyzed in detail. 
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Fig. 2. The refrigeration heat ratio r R versus the temperature ratio of cycle curves for the possible surface area ratios at given r v = 3 and X = 0.16; 
(a) 1 < r A < r v , (b) r v < r A < r v [l + 2y/nV L (l + y/z) In r v /l c (T H )A L ]. 



Fig. 3. The refrigeration heat ratio r R versus the volume ratio of cycle curves for the possible surface area ratios at given z = 0.8 and X = 0.16; 
(a) 1 < r A < r v , (b) r v < r A < r v [l + 2 y/nV L (l + y/r) In r v /L c (T H )A L ]. 


r £ 


In following paper, for convenience of analysis, one may introduce the refrigeration heat ratio r R = Rqb/Rc, the COP ratio 
= Sqb/sc and the dimensionless thermal wavelength k = • For a given k = 0.16 and some different surface area 


ratios in the two ranges of 1 < r A < r v and r v < r A < r v [l + 2y/jrV L (l + y/r) In r v /L c ( T H )A L ] in micro/nano scaled 
system, the refrigeration heat ratio r R versus the temperature ratio (at given r v = 3) and volume ratio (at given r = 0.8) 
of cycle curves can be plotted, respectively, as shown in Figs. 2 and 3. It is clearly seen from Fig. 2 (a) in the case that the 
refrigeration heat ratio r R increases as r increases at given other parameters. Especially, the value of r R is larger than one 
and as a decreasing function of the surface area ratio r A as r > 0.25 but smaller than one and as an increasing function of 
the surface area ratio r A as r < 0.25. Moreover, for a given value of 1 < r A < r v and r = 0.8, one can observe from Fig. 3(a) 
that r R is not a monotonous function of r v and exists a maximum value at the special value of r v . 

We can find from Fig. 2(b) and 3(b) that the refrigeration heat ratio r R is less than one in the case of r v < r A < 
r v [l + 2y/rzV L (l + yfr) In r v /L c ( T H )A L ] and increases as r and r v increase at given other parameters. Further, in the 
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Fig. 4. The coefficient of performance ratio r £ versus the temperature ratio of cycle curves for the possible surface area ratios at given r v = 3 and k = 0.16; 
(a) 1 <r A < r v , (b) r v < r A < r v [l + 2 yfnV L (l + y/x) In r v /L c (T H )A L ]. 



r v 

Fig. 5. The coefficient of performance ratio r e versus the volume ratio of cycle curves for the possible surface area ratios at given x = 0.8 and k = 0.16; 
(a) 1 <r A < r v , (b) r v < r A < r v [l + 2 *JnV L (l + ypt) In r v /L c ( T H )A L ]. 


case the refrigeration heat always decreases as the surface area ratio r A increases. Figs. 4 and 5 give the curves of the 
COP ratio r s of cycle varying with the temperature ratio (at given r v = 3) and volume ratio (at given r = 0.8) of cycle 
for different surface area ratio in the range of 1 < r A < r v [l + 2yfnV L (l + ^/r) In r v /L c (T H )A L ]. The COP ratio r £ of 
cycle is a monotonous increasing function of the temperature ratio r and the surface area ratio r A but not a monotonous 
function of r v , and exists a minimum value at the special value of r v at a given range 1 < r A < r v . For the case of 
r v < r A < r v [l + 2y/jtV L (l + y/x) In r v /L c (T H )A L , the COP ratio r £ increases as r or r v increases but decreases as r A 
increases at given other parameters. In particular, the COP ratios r £ of cycle in two cases are always less than one and the 
COP of the cycle may be improved when the volume ratio r v approaches at the value of the surface area ratio r A then the 
COP of cycle equals the ideal Carnot COP s c . It is noticed that the results obtained here based on the boundary effect of gas 
particles come about when the de Broglie’s wavelength is comparable to the dimensions of the system. 
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4. Conclusions 

In the paper, by considering the quantum boundary effect of gas particles, which appears if the system is confined in 
a finite domain, the performance characteristics of a Stirling refrigeration cycle in micro/nano scale are reanalyzed. It is 
found that a micro/nano scaled Stirling refrigeration cycle working with an ideal gas does not possess the condition of 
perfect regeneration due to the existence of quantum boundary effect in a finite domain. Furthermore, considering the 
boundary effects, the refrigeration heat and COP of the micro/nano scaled Stirling refrigeration cycle are derived. They are 
the functions of the temperature ratio r of the two heat reservoirs, the volume ratio r v and surface area ratio r A of two 
isochoric processes and so on. In particular, the influence of the boundary of domain i.e. surface area on the performance of 
a micro/nano scaled Stirling refrigeration cycle is analyzed by introducing r R and r £ . It is found that the refrigeration heat 
of the cycle may be large than the refrigeration heat of the macro scaled Stirling refrigeration cycle as 1 < r A < r v . That 
is, more heat is absorbed per cycle from cold reservoir. In spite of the fact, the COP of the cycle cannot exceed the COP of 
the macro scaled Stirling refrigeration cycle s c , since the corresponding work input per cycle must increase. However, as 
mentioned in Ref. [3], the number of cycles for refrigerator is an important performance parameter. In a real refrigerator, 
lost work also increases and COP decreases with an increasing number of cycles due to the friction losses. Consequently, it 
is understood that in the case of 1 < r A < r v the quantum boundary effect can be advantageous for a micro/nano scaled 
Stirling refrigeration cycle since the same refrigeration heat is obtained with a lower number of cycles. In the opposite 
direction i.e. r v < r A < r v [l + 2jnV L (l + ^/r) In r v /L c (T H ) A L ], the refrigeration heat of the cycle is always smaller than 
the refrigeration heat of the macro scaled Stirling refrigeration cycle. Further, in the case the corresponding work input 
per cycle decreases and the COP of the cycle is also always smaller than the COP of the macro scaled cycle s c . Thus, in the 
case, the quantum boundary effect is always disadvantageous because it causes not only an increment in the number of 
cycles for given refrigeration heat but also the decrement of the COP of the cycle. Of course, for the present micro/nano 
scaled cycle, the COP of the cycle may be as large as possible and approach the ideal Carnot COP s c , by selecting properly 
the surface area ratio r A of cycle approaching the volume ratio r v of cycle. The results of the paper are very important 
and useful for designing of the micro/nano scaled Stirling refrigeration cycle. In addition, when selecting the values of 
surface area ratio r A of cycle deviating from the value of volume ratio r v of the cycle, the COP of cycle decreases largely. 
Thus this case leads to confirming experimentally the quantum boundary effect in micro/nano scaled devices. 
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